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Abstract: We consider the equilibrium partition function of an ideal gas of Dirac 

fermions minimally coupled to torsion in 2 + 1 dimensions. We show that the energy- 
momentum tensor reproduces the Hall viscosity and other parity violating terms of first 
order in the torsion. We also consider the modifications of the constitutive relations, and 
classify the corresponding susceptibilities. An entropy current consistent with zero pro¬ 
duction of entropy in equilibrium is constructed. 
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1 Introduction 

The description of transport phenomena and dynamical response in the framework of ef¬ 
fective field theory has experienced great progress in the last few years. This is due, to 
a large extent, to a better understanding of the role played by the quantum anomalies 
which underlie new macroscopic parity-violating effects at low energy [1], In this respect, 
the construction of the thermal partition function as a derivative expansion of a time- 
independent background has turned out to be an important tool to get information about 
the non-dissipative part of the constitutive relations of hydrodynamics, without using an 
entropy approach [2, 3]. In relativistic systems the first order of such expansion usually 
has a relatively simple structure connected with anomalies [4-6]. 

For non-relativistic systems there is also considerable interest in establishing the precise 
connection between the partition function and Hall transport [7-12], While in these systems 
the use of torsional Newton-Cartan geometry appears as natural, the role of torsion in a 
relativistic setting is less clear [13, 14]. Originally the torsion-dependent effective action of 
massive Dirac fermions at zero temperature and density has been examined in great detail 
in refs. [15, 16] with the focus on the renormalization effects on the Hall viscosity. Here we 
are interested in the application of the methods of [2, 3] to the thermal partition function 
to linear order in the torsion. 

The main purpose of this paper is to compute, to linear order in the torsion, the stress 
tensor and charge current that follow from the equilibrium partition function for a Dirac 
field minimally coupled to torsion. A first result of our computation is that the relationship 
between the Hall viscosity and the spin density, f/ = (£)/2 = — ('F\E f )/4, naturally appears 
as an equilibrium susceptibility relating the stress tensor with certain parity-violating 
combination of components of the torsion tensor. Other results are the explicit construction 
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of an entropy current consistent with zero production of entropy at equilibrium, and the 
derivation of new susceptibilities relating the charge current and the pressure with torsion- 
dependent vector and scalar data. 

The organization of the paper is as follows. In section 2 we present some details about 
the notation, and the partition function we use in the subsequent computations. The form 
of the stress tensor at thermal equilibrium is given in section 3 in terms of combinations 
of scalar, vector, and tensor background data depending on torsion. The modification in 
the Landau frame of the constitutive relations in the presence of torsion is presented in 
section 4, as well as an expression for the entropy current compatible with zero entropy 
production. We conclude in section 5 with a short summary of our results. 


2 Partition function of a Dirac fermion coupled to torsion 


We begin with the action for a Dirac field, 


S= d 3 x det e“ 


+ m'F'F 


7^0z) = e^(x)7 a , (2.1) 


in the most general static background with torsion 

ds 2 = G^ u dx^dx” = — e 2cr ( x ) (dt + at(x) dx*) 2 + g l] {x)dx l dx 3 . 

A = Aodt + Aidx 1 1 

T a = -T^dx' 1 A dx u . 


This geometry with torsion has been also considered in other studies of the torsion re¬ 
sponse [17]. Note that we are using the prescription of minimal coupling, where the action 
depends on the torsion only through the covariant derivative. The notation of ref. [18] has 
been adopted, where the Dirac adjoint is defined as T = Tfry 0 . For the Dirac matrices 
7“ we choose the representation {7°,7 1 ,7 2 } = {—io 3 , 02, — 07}. The torsion 2-forms T a 
are specified in terms of the frame field e“ and the torsion tensor by T^ iu a = T /Jiy CT e“. The 
covariant derivatives 


^p. ^ = (<9p + juj^ ab 7 ab — L4 At ) 1 F, 7 ab = \ [7 a, 7b] 


(2.3) 


are written in terms of the spin connection 0 J^ b , which is specified by the 1-forms u ab = 
u^dx^ appearing in the Cartan structure equation 

de a +ix a b Ae b = T a . (2.4) 


By defining the contortion tensor, antisymmetric in the last two indices, as 


(2.5) 
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it turns out that the spin connection can be expressed in terms of the Christoffel symbol 
T° (g) and the contortion tensor as follows 


LU 


ab 


= —e^d^e® + e bv - K^) e“ 

= ^ ab (e) + I<^ 


ab 


( 2 . 6 ) 


(The square brackets denote antisymmetrization of indices Ar fW i = ^{A pL , — A vp ).) This 
formula contains the affine connection T" w = — K^J 7 in the presence of torsion, 

which has an antisymmetric part given by 


p<T _ p(T _ rp O 
1 11V X Ufl flU 


(2.7) 


The part of the spin connection in the absence of torsion which is called the Levi- 

Civita connection, is uniquely determined by the frame field. It is obtained from (2.6) by 
setting K^J 7 = 0. 

In order to vary the metric and torsion variables independently we will assume, ac¬ 
cording with ref. [19], that the torsion tensor with the last coordinate upper index T a 
and the metric components G are independent of each other. This will be important 
soon. Note that in ref. [14] the contortion is considered as an independent variable, instead 
of the torsion. 

With the previous relations and 


{7c Tab} = -2e cafe , e 012 = 1 =-e 012 , 

e^ = dete d a G bc e»y h e^ 


( 2 . 8 ) 


we can see the effect of the torsion on the three-dimensional Dirac action. The substitution 
of (2.6) and (2.7) into (2.1) reveals that S may be viewed as the action of a Dirac field in 
a torsion-free background, but with a modified mass. The new mass rh, that now depends 
on the torsion, is obtained by the replacement 


m —>■ rh = m + 5m = m — ^ flup T flu x G\ p . 


(2.9) 


It is remarkable that a generic non-minimal coupling of the fermion to torsion may be 
included simply by replacing 5m by £5m, where £ is a free parameter [16]. 

In this way, the partition function at zeroth derivative order in the variables {A p , G llUl m) 
is written as 


/' 


/ d 2 x^—V(T,^rh), 


T = T 0 e~ 


g = A 0 e~ 


( 2 . 10 ) 


where T 0 1 is the period of the imaginary time, and the function V is the pressure in terms 
of the temperature and chemical potential. Therefore, we expect that the functional 

1 


W[G pu ,T x ] = - 


d 2 x^-(W)e^T x G Xp , 


( 2 . 11 ) 


contains all the information about the static linear response to lowest order in the derivative 
expansion of the torsion. Here we have used the fact that the spin density is related to 
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the pressure by dV/dm = (\F'F), which follows from the form of the Dirac Lagrangian, 
eq. (2.1). 

Although no specific form of V is required in the following computations, we write 
down the expressions of the pressure and the spin density at non-zero temperature and 
density in the free-held case, 


2 vr 

Li 3 

\m\ 

T 2 

2ir 

mT 

In 


H — \m\ 


Li 2 — exp 


T 

fj, — \m 


+ Li 3 - exp 


-/r — \m\ 


2ir 


In 1 + exp 


T 

/r — \m\ 


+ Li 2 exp 
+ In ( 1 + exp 


~H — \m\ 
T 

H — \m\ 


( 2 . 12 ) 

(2.13) 


where Li n (x) is the polylogarithm function. 

3 The equilibrium stress tensor at linear order in the torsion 

We can now find the energy-momentum tensor 0^ which follows from the partition func¬ 
tion by differentiation with respect to a, aj and gd , with the understanding that is 
independent of the metric. From eq. (2.11), the variational formula for 0 qo leads to 


@oo = - 


= —e 


T 0 e a 5W 
y/9 2 dv 
2ct ( dj**) 

V da 


+ (TT) 5m - 


d5m 

~d^' 


(3.1) 


The most obvious way to obtain <9('F'F)/<9 ct is to consider the energy density rewritten in 
the form 


^ BV BV 

£ -~ v + T W' + ^ 

da' 

Then differentiating with respect to m yields the relationship 


(3.2) 


= _/**> * a < J,,1 '> 
dm K ' da ' 


(3.3) 


and eq. (3.1) becomes 




Ola 


Be 


= e““ + 0^) ^ y ^(Toij - cijToio) - e Aa 5m 


(3.4) 


In the computation of the last derivative eq. (2.9) has been used, while keeping T ( tI/ A 
fixed. Here e 12 = 1 /y/gi- Note the invariance of this result under time reparametrization, 
t —>• t+(j)(x), x —>• x. In the light of (3.4), it is convenient to introduce another pseudo-scalar 


- 4 - 




















quantity (besides 5m) constructed from the torsion as E = — u p ' T lwp e vpx u x . This quantity, 
evaluated for the equilibrium fluid velocity, u^ = <5Qe _cr , becomes E = —e _c7 e u To ij, where 
To ij = Toij — ajToio. Thus the change of ©oo induced by the torsion reads 

@00 = e 2<7 —— 5m + e 2a { \E'\E f ) ( — - 5m 
dm \ 2 

The other components of the energy-momentum tensor are easily computed: 



_ Toe"' 7 5W 

s/ 92 5ai 

= <**> (x £ ‘ J ( r J - «j r o) - • 


(3.6) 


where T ^ is the torsion vector defined by T p = T /. This expression suggests the intro¬ 
duction of two pseudo-vectors orthogonal to u p given by 


X p ' = -e^ p u„u x T Xpa u a , 
m = -e pup u u T p . 


In equilibrium these quantities evaluate to 


X* = e- 2a e^T 0j0 , 

W { = e * J (Tj - OjTq). 


(3.8) 


For the equilibrium stress tensor the functional differentiation of W yields 


Qij 


2T ° e g a ik a jm SW 

y/92 5g km 

-^(^)(e ik ‘X ok mg mj +e j Ho km g mi ). 


(3.9) 


Remarkably, this expression may be expressed in terms of the pseudo-tensor given by 




c vX Pi 


• )•> 


(3.10) 


where is the shear tensor, 


= A pa A ^(V a ^ + V 0 u a - G a pV p u p ), A pa = G pa + u p u a . (3.11) 


While in the absence of torsion the equilibrium values of a pu and a pu vanish, it turns out 
that when torsion is present, the covariant derivative of the equilibrium velocity u p = 5 q e~ a 
acquires a non- zero contribution proportional to the torsion, which now is included in the 
affine connection. Correspondingly, the shear tensor at equilibrium is given by 


— 


(Kap X + Kpa ) M A + G a pK uX u U 
—(T Xa p + T X /3 a )u x + G a pT Xv v u x 


U x = 5^e~° 


(3.12) 
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Hence, apart from a 13 no more pseudo-tensors are needed to express the equilibrium value 
of ©*U Thus eq. (3.9) may be rewritten in the form 


0 b' = I(^) (a ij -Eg ij ) 


eq 


(3.13) 


In a general configuration eq. (3.12) is not longer satisfied, and in addition to a pu , it 
is convenient to define another pseudo-tensor related to torsion given by 


TflV 

r (2) 




(3.14) 


where t p ^ = — A pa A v P{T\ a p + T\p a — G a pT\)u x . This will be important when we consider 
a possible generalization to non minimal coupling in the last section, because then will 
appear in the constitutive relations. 

It is also possible to derive these results directly from the energy-momentum tensor 
obtained from the action (2.1), without using the partition function (2.11). With the 
assumption that the spin connection does not depend on the frame field, the functional 
derivative of the action with respect to e p produces the non-symmetric tensor 

T ,pv = -- e ^g pv = - V V)^- (3-15) 

(3 O&a A 

The symmetry may be restored by adding the tensor AT pv given by 

A T pv = ^(V A + T x ) (T{ 7 A/ \ 7 "}T), (3.16) 


where the covariant derivative is performed with the affine connection including the con¬ 
tortion, and T\ denotes the torsion vector T\ = T pp p . Using the equations of motion this 
gives the usual form of the stress tensor 

Q pu = T niv + AT ^v = + (^ o i/))T. (3.17) 

To obtain the dependence of (Q pu ) with the torsion, we note that the contortion tensor 
in the covariant derivatives always multiplies the term (TT). A second kind of contribu¬ 
tion comes from the fact that the fermion field satisfies a Dirac equation with a torsion- 
dependent mass m. As a consequence the corresponding perfect fluid constitutive relation 


(&"') = e(T, fi, m)u p u u + V{T, /q m)A py , (3.18) 

contributes to the expectation values of ©oo and Q lJ with terms proportional to de/dm 
and dV/dm respectively. Then the explicit evaluation of K^ ab in the spin connection shows 
that the equilibrium values may be written as the combinations 

(©oo) = (**> - e 2a Sm)j + e 2 CT |^m, (3.19) 

( 0 b) = _(* *} f^^ okmg mj + e 3k % 0km g mi ) + + ^5mg i3 . ( 3 . 20 ) 

These are in perfect agreement with (3.4) and (3.9). The result for (0q) coincides with (3.6), 
and does not include contributions of the second type because they are not generated by 
the perfect fluid constitutive relation. 
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3.1 Kubo formula for Hall viscosity from a contact term 

So far we have been working at the level of the partition function, but it is instructive 
to consider the response to a general time-dependent applied torsion in the framework of 
linear response theory. With this in view, we now consider eq. (3.20) as the static limit 
of a Kubo formula for the stress tensor response to an externally applied time-dependent 
torsion (not necessarily uniform). In order to obtain such a formula, we assume that the 
only non-zero components of the torsion tensor are Tnj k = —Tj ofc, and we also consider 
that the metric is Minkowskian. Thus the change of the Hamiltonian to linear order in the 
torsion becomes 


H\ = — j d 2 x 1 F'I'<5 m 


= -- d 2 xWe l T 0ij , 


(3.21) 


The induced change of the expectation value of 0jj is obtained from linear response theory 
in the form 


S(@ij(x,t)) 


d 2 


/ 5eij(x,t) 
V \fiT okm (y,t) 


0 


+ 7 / dfe-tf*-*) 
4 ./__ 


d 2 y TT(r/, t )]) 0 e km T okm {y, t), 


(3.22) 


where 77 —>• 0 - * - , and the subindex 0 means an equilibrium average with respect to the unper¬ 
turbed Hamiltonian. The first term (contact term) has its origin in the explicit dependence 
of 0 with the torsion, and corresponds exactly to the first summand of eq. (3.20), 

/ SQij(x, t) \ = 1 /^ o / e ik 5 mj + e jk 5 mi + e km S ij\ _ y ) _ (3 23 ) 

\ cTq / 0 ^ ' ' 

The second term involves the retarded stress tensor-spin density correlator defined by 


3 ^\j (*^7 


- l - lim 0(t) (J@ij(x, t), TT(0, 0 )]> o e 

Z 77—^O -1- 


(3.24) 


We can express eq. (3.22) in terms of a response function to the torsion Xij km{q, w), which in 
the Fourier domain relates d(0ij) = Xij kmTokm- Therefore eq. (3.22) leads to the following 
Kubo-type formula 

Xijkm(q^) = -\(mo (e ik 6™i+ei k 6 mi ) - + i(TT) 0 ^)e fcm , (3.25) 

where 

POO P 

Yij(q, ui) = J dte^+^J d?x e~ iq x ([0ij(*, t), TT(0,0)]) Q . (3.26) 

The response function Yij(q,oj ) is similar to the (integrated) stress-strain form of the re¬ 
sponse function in ref. [ 20 ], where the role of the spin density is played by the antisymmetric 
part of the strain generators J af 3 . 
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We see that the Kubo formula (3.25) is consistent with the static result of eq. (3.20) 
only if the ordered limit lim^o hm^-^o Ty(g, oj) satisfies 

e W 

lim lim Yij(q,u ) = lim Yij(q,uj = 0) = —^ —. (3.27) 

g—>0 oj—> 0 g—>0 2 (7m 

This requirement can be expressed as a thermodynamic sum rule valid for q ^ 0, 


lim 

q—> o 



doj lmYij(q , cu) 

7T w — i0+ 


V dV_ 

2 9m’ 


(3.28) 


which is similar to other sum rules giving different susceptibilities. It follows that the static 
response to T^m is determined by a part of the contact term, namely 

lim Xijkm(q, w = 0) = f, U k 8 m i + fJ k 5 m ) , (3.29) 

q—> 0 \ / 


so that 

fj = - lim SijekmXij km(q, w = 0). (3.30) 

4 g—►() 

The coefficient fj = — (T'I')o/4 will be identified below with the Hall viscosity. 


4 Torsion-dependent constitutive relations 


We now consider the corrections to the constitutive relations imposed by the partition 
function (2.11). The argument essentially follows ref. [2], As these corrections involve 
the expectation value of the particle current = ilLy^T, we also need the modification 
induced by the torsion on the equilibrium current. This is simply given by 


T 0 e a 5W 
y/92 SA 0 


-e a ^-5m, 

am 


= 0 . 


(4.1) 


This may be immediately obtained from the relation between the particle density and the 
pressure, n = e a dV/8 Aq, which implies 


4.1 


9(TT) 

Stress tensor and charge current 


dn 

dm 


(4.2) 


The most general parity violating modification of the currents that includes torsion must 
be written in terms of the non-zero quantities in equilibrium 5m, S, W and fj ^ v . Then 
the non-dissipative parts of the constitutive relations in the Landau frame take the form 


T» v = £u^u u + VAT - (xrSm + x~E) A^ -fja^, 
J' 1 = nvf 1 + xwW 11 + xxW, 


which involve four new transport coefficients xr, Xs, XWiXx and the Hall viscosity fj. 1 To 
find them, we specialize to equilibrium by using u^ = <5ge _cr to zero order in the torsion, 

1 The term (xrSm + xeE) A pl/ would be replaced by (; \bB + xnS! + xr 8m + XbH) Ain the complete 
constitutive relation, where B = 1 e p ' 1/f> F l/p and the vorticity is defined by fi = e MI ' p M M chu p rather than 

e pvp u ll V„u p . Then the effect of the additional term related to the antisymmetric part of the affine connection 
in the presence of torsion, AS! = —e p " p u p Tl P H A , is considered separately in the coefficients xr and xs- 






and Aoo = &oo = Xq = Wo = 0. Because the torsion-induced corrections can be expressed 
in terms of changes in the fluid variables (8T,5fj,,5u l ) according to 


STq = -e a (e + V)Su i , 

5f = n5u i + xwW* + Xx X\ 

f dV dV \ 

5T lJ = <5T + ~ XT 8m - xsSj g v - f) a lJ , 


(4.4) 


we can equate these first three equations to (3.5), (4.1) and (3.6). This gives 


de dn 

dfdg 
de dn 
dfdg 


de dn 
~dJi~dT 


5T = 

dn / 

) = 

dn \ 
dn I 

dr \ 


' de T .\ de dn ] (T4 ’) dn„ 

O” 7d/i0mJ m 2 9/1 1=1 ’ 

f de T .\ de dn] <9n„ 

~d^ + ( ) ) + dfd^l 6m + ^~dT : 


5u l = 



e + V 


(4.5) 


We now determine the transport coefficients appearing in the constitutive relations by 
equating the last two equations in (4.4) with (4.1) and (3.13). With eqs. (4.5) and the 
thermodynamical derivatives 


dV / dV dn dV dn\ f / de dn de dn \ 

~de = \dT~dn ~~ 9 / 7 dr) / \dfdjl ~ ~dndTj ’ 
dV _ / dV de dV de\ / f de dn de dn \ 
~dn = \Wdfi + 9 ) 7 df) / \dT~dJl ~ d^dTJ ’ 


we are left with 


« = < 4 *> (‘ - °i ) • 

xh = -W (~ + \ 


V = 


(TT) 


1 /=F T\ n 
1 / T T\ 71 

Xx = - 2 (^) 7T p- 


(4.6) 


(4.7) 


These expressions, which uniquely determine the new transport coefficients in term of the 
angular momentum density, are one of the main results of this paper. 
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4.2 The entropy current 

Finally we consider the contribution to the entropy current which arises from the partition 
function in eq. (2.11), 

W = j d i x ^ 2 C ^ dV{T d ^ m) 5m, dm = -\e^ p T, up . (4.8) 

If s denotes the entropy density that follows from the thermodynamic potential at zero 
derivative order, s = dV/d T, the standard formula for the entropy yields the contribution 


= d(T 0 W) 
dT 0 



d 2 V . 
dm dT Sm 



(4.9) 


Therefore a condition that must be fulfilled by the correction to the entropy current dJg is 



(4.10) 


Let us consider the following tentative expression for the entropy current containing the 
dependence on the torsion: 

= su 11 - ^ (xwW p + xxX' 1 ) + (6= H + b r dm)u p + d w W p + d x X p . (4.11) 

It follows then from this form of the entropy current that the four coefficients b =, b~ f, dw 
and d x are determined, in a unique way, by the condition 


dJ$ = sdu° + e~ a 


ds XT , ds A 


T 



f) c 

= e~ a ^-dm. 
am 


+ (fe 


+ bj;dm?)u p + d\\rW p + d x X p 


Then one finds that 


6= 




dw 


d x 


(4'4'} 

2 T ’ 
2(4/'I') 
T ’ 
(4-4') 
4T ’ 
(4/4/) 
2T~ 


(4.12) 


(4.13) 


Note that in order to obtain these results it is necessary to use the values of dT and 
dfi given by eqs. (4.5), as well as the expressions for du° = — a^du *, IT 0 , X°, and the 
thermodynamical identities 


ds 1 
de = T’ 


ds /j, 
dn ~ T 1 


e + V = Ts + jin. 


(4.14) 
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Thus one has 


(4.15) 


ds 1 de n dn 
dm T dm T dm 

It is remarkable that the results in eqs. (4.13) imply that the spatial part of 5Jg is 
identically zero at equilibrium: 


<54 = aSu‘ + (-£* *«• + iw) W‘ + (~XX + ix ) X' 
e + V l 4 


(4.16) 


= 0 . 


Thus the condition of no entropy production is satisfied since the divergence of 5Jg is 
necessarily zero under these conditions. 2 


5 Conclusion 

In this paper, we have studied the modifications of the stress tensor of an ideal gas of Dirac 
fermions coupled minimally to torsion. We have concentrated on the part coming from 
the equilibrium thermal partition function to linear order in the torsion, which in 2 + 1 
dimensions may be easily expressed in terms of an effective mass shift. Differentiating 
with respect to the metric we have obtained different types of background data linear in 
the torsion, including two pseudo-scalars, two pseudo-vectors and a pseudo-tensor field. 
We have also found an expression for the entropy current consistent with zero entropy 
production at equilibrium. It turns out that the non-dissipative part of the constitutive 
relations involve five susceptibilities, four of them expressing the modifications of the stress 
tensor and the current in response to the applied torsion. Because of the way the torsion 
affects the covariant derivative of the fluid velocity, the remaining susceptibility may be 
written as the proportionality factor that expresses the equilibrium response to the parity 
odd shear tensor. Therefore it would be identified as the Hall viscosity. 

Throughout this paper we have assumed that the fermions are minimally coupled to 
torsion. Considering non-minimal couplings introduces an additional parameter on which 
some of response functions may depend. Specifically, one has to make the replacement 
5m —>• £5m for arbitrary £, therefore leading at first sight to fj —>■ £fj . This has been used 
to question the existence of a connection between Hall viscosity and torsion response [14]. 
However, the analysis sketched above ignores the fact, pointed out in section 3, that there 
are actually two independent tensors, given by eqs. (3.10) and (3.14), which coincide at 
equilibrium (see eq. (3.12)). Therefore, both tensors must be used in the modified consti¬ 
tutive relation 

= -fj cr^v - m t 1 ^ • (5.1) 

2 If torsion is present, the covariant differentiation is defined with respect to the non-symmetric connection 
rJD(s) — i/ 7 . However, it is possible to define a modified divergence (see e.g., ref. [19]) according to 

which for a vector field W reduces to the torsion-free result = {\/—g) 1 d—gV 11 ). 

* 

Thus the conservation or the entropy current found above must correspond to Jg = 0. 
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Now, as the covariant derivative in the presence of torsion involves the affine connection 
T° = T° u (g) — K jUJ a , rather than r a ^ v {g) — , the parity odd shear tensor in the 

last term of eq. (4.3) remains unchanged even if £ ^ 1, and the same applies to rj which 
cannot depend on £. Therefore the torsional response in the presence of non-minimal 
coupling requires that 

V2 = (l-0^- (5.2) 

In this way the Hall viscosity remains intact, and in the absence of torsion is the only 
non-dissipative term in the constitutive relation. 

Finally, it would be desirable to combine the results from holographic models with 
those presented here, but it does not appear to be an easy task. The main reason is 
that there is no unique way to produce Hall viscosity and angular momentum density in 
holography [4, 21-23]. For example, the holographic model considered in ref. [4] gives rise 
to nonzero angular momentum density but vanishing Hall viscosity, while for the model 
considered in ref. [22] the ratio between these quantities is compatible with the universal 
value 1/2 from field theory. In this respect, it may be interesting to note that for a gapped 
system, such as a free massive Dirac fermion, the spin density of eq. (2.13) is nonzero 
at n = 0, while the angular momentum density (£) from certain holographic models [4] 
vanishes at fi = (1. Thus the relation (£) = \dV/dm that determines many features of the 
results in this paper does not seem to be satisfied in those models. Further studies are 
required to understand all these points. 
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